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Abstract 

We calculate distributions for r+r^7 production at LEP 1 taking into account a potentially 
existing anomalous magnetic moment of the r lepton. The existing upper limits for |a^| 
are known from the dependence of the decay t^t^'-/ on and are of the order of 

(1 — 5)%. We show that such limits are also sensitive to linear terms in a,-, which are of equal 
importance at |a^| ~ (1 — 2)% and dominate below this value. Contributions from an electric 
dipole moment d^- do not interfere with the electromagnetic vertex or with the anomalous 
magnetic moment. Appropriate formulae are derived. 



1 Introduction 



There are few basic physical parameters which may be measured best at LEP 1. One of them 
is the anomalous magnetic dipole moment of the r lepton. 

The matrix element of the electromagnetic current of a particle with mass m and spin | 
has the form illl: 



eu{p' 



uip), (1) 



2m 

with q = p' — p and cr^'^ = | [7^^, 7^^]. The corresponding magnetic dipole moment is given by: 

f^ = ^[Fm+F2m, (2) 

2m 

with -Fi(O) = 1. The anomalous magnetic dipole moment of a fermion / is defined as 

af = FU0) = ^f^f - 1. (3) 

Without radiative corrections or anomalous contributions, it is aj = 0. The lowest-order 
QED correction is af = a/(27r) = 0.001161 0. The magnetic moments of electron and muon 
are measured with high precision with the spin-precession method 0]: 

eh 

fie = (1.001159652193 ±0.000000000010) , (4) 

2me 

/i. = (1.001165923 ±0.000000008)-^. (5) 

2m^ 

They are in complete agreement with higher-order theoretical predictions and represent an 
important test of the Standard Model 0. 

For the short-lived r lepton the spin-precession method cannot be used. To our knowledge, 
the first measurement of the form factor -^2(1?^) of the r lepton was performed with the e^e~ 
annihilation experiments at PETRA 0. The result was ^2(9^ 7^ 0) < 0.02 (95% CL) for 
values of q^ in the range (5 — 37 GeV)^. The best limit for F2{q^) is due to an analysis of 
Escribano and Masso [§]. Using an effective lagrangian approach to the Z° width, they derived 
from experimental values for the frZ° couplings the limit -F2(M|) < 0.005 (95% CL). This is 
one order of magnitude larger than the theoretical prediction 

a*^ = 0.001177 (6) 

However, since the g^-dependence of F2 is unknown, it is impossible to infer the static magnetic 
moment from the above measurements of -^2(0'^) at q^ 7^ 0. 

The first direct determination of the anomalous magnetic moment of the r lepton, i.e. of 
^2(0) is due to Grifols and Mendez fig (95% CL): 

|a^|<0.11 (7) 

They derived this limit from a calculation of the quantitative effect of the electromagnetic 
current (0) for on-mass-shell photons (with q^ = 0) on the partial decay width: 
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10247rm^ sin 6w cos^ 9w 
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The standard model contribution is Fq and v = 1 — 4sin^ ^^w^^ a = 1. 

The hmit (|^) has been improved recently by a dedicated analysis of the L3 collaboration [|11 
[g (95% CL): 

\ar\ < 0.049 (9) 

L3 analyzed the distribution 



aE^ IMirMzm^ sm Ow cos"' Oi^ ' i" j 

As may be seen from the above expressions, both determinations of a,- neglect terms linear 
in a,-. A simple estimate based on the square of matrix elements indicates that the linear 
terms are in fact suppressed in comparison to terms with by a factor r(s) = m'^/s, with 
s being the scale of the process. At LEP 1, s = M| and r(M|) ^ 1/2633 = 0.38 per mil. 
Thus, as long as a limit on \ar\ considerably exceeds this value, the neglect of linear terms in 
ttr seems to be justified. 

Recently, a calculation of distributions for the process 

e+e" ^ (Z°,7) ^ r+r^7 (11) 

at LEP 1 energies gave hints on a stronger influence of the linear terms in than the above 
estimate based on r(M|) suggests |T^. Because the argumentation is based on a numerical 
calculation, one may ask for an analytical re-calculation in order to better understand whether 
the effect really exists and what its origin is. Last but not least one needs an estimate for the 
range of values of a,- where an improved calculation is necessary. In addition, besides the total 
rate and the energy spectrum of decaying bosons, the angular distribution of photons and 
T leptons in the center-of-mass system (cms) may be well measured at LEP. Thus, one should 



clarify whether these distributions may be of use for a determination of limits for a^- [0 



In this article, we calculate several distributions and the total rate for the reaction ([Tl|) 
arising from the Feynman diagrams of Figure ^ with account of the anomalous part of ver- 
tex ([^). We formulate the problem exact in m-r, but will later neglect higher powers of rrir 
where they are evidently superfluously for numerical estimates. Further, we can safely ne- 
glect anomalous contributions arising from the initial-final state interferences, from the 
interference, and from 7 exchange. Other Standard Model contributions like e.g. initial-state 
radiation are assumed to be known. 

In Section H, we discuss the analytical formulae for d^(T/dcos6'^di?^, da/dEy, atot ■ Mod- 
ifications arising from an electric dipole moment of the r lepton are mentioned. Section ^ 
contains a numerical discussion. The application of simple cuts is also considered. In Ap- 
pendix A the phase space parameterization and the choice of four-vectors are explained. 
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2 Distributions for e^e 



Z 







We used form for the calculation of the square of matrix elements, then inserted the 
products of four-momenta in the center-of-mass system as determined in Appendix 0, and 
finally performed the various integrations. The twofold cross-section differential in energy and 
production angle of the photon in the cms is: 



dV 



da:'d cos 9^ 



(.^ + a^)llni±|(l + cos^^, 
a2[(l-x')^(l-cos2 e^) + 8x'] 



m. 



2«r 



+ 



3{v^ + a^) 
8(1 -x') 



:i + cos' ^^)(l + x'')-ln- 



P' 



(1 + cos^ e^){i-x'y -8x' cos^e. 



4a' 



121 



c 



aGlMl 



vr 967r |s - M| + iM^r^(s) |^ 



Ami 
x's 



and x' is the normalized invariant mass s' of the r+r : 



X 



s 
s 



It is related to the center-of-mass photon energy E-^: 



l-x') 



(13) 
(14) 

(15) 
(16) 



In (|12|), some terms proportional to a^- are enlarged compared to due to a logarithmic factor 

In ~ ln(x's/m2), which is 7.88 for x' ^ 1 at = Mz = 91.187 GeV and = 1.777 

GeV. The simple numerical factors also disfavor compared to a,-. Thus, the suppression 
factor r(M|) overestimates by far the relation between the quadratic and linear terms. We 
should also remark here that the logarithmic term ■ ln(s/m^) does not represent a mass 
singularity. It arises from the combination (^i • ar/mr) ■ (m^/^i) ln(s/m^), where ar/rrir 
has its origin in (^ and the factor {mr/^/s)\n{s/m'l)) arises from photonic corrections and 
vanishes for small rrir/ ^/i as it should do. 

As one would expect, distribution (|12D is smooth in cos^ and independent of cos^^. 
The Born angular distribution of the r lepton is smooth and the final state photons are 
preferentially coUinear. The forward-backward asymmetry (in cos^-y) must vanish and thus 
odd powers in cos 6'^ are forbidden. 
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The integration of ([T^) over cos 9^ is simple and yields 



dcT 
dx' 



P'Cl (1 - x') 



+ 4(1-^') 



l + (3' 



(3' 1-/3' 



7 -3a' 



'1 — X 



r\2 



+ {v^ + a^)x' 



m. 



2^r 



+ 



1 -x' ' 



1 + X' 



'2\ 



In 



1 + 



- 1 



1-/3' 

The terms proportional to may be found in (|l^) already. 
An integration of ([1^) with lower cut x on the photon energy, 



(17) 



<^hard(X) 



dx'^^. 
ax' ' 



with 



yields 



X 



1 

+ 24 



{v^ + a^)\ix-2)\n^ + 2--] - 3a\x - 2) 
\ 2 



(19) 



^2^_0_J^ + (v^ + a^)x'{3 - 2x) 



s 2 



2 T 



+ iv' + a') 



+ 21n 



-x^ - -(2 + x) In ^ + 2x + 2 Li2(l 
4 2 m: 



X 



2Li2(l) 



1 — X 



In^ - 1 



(20) 



The function Li2 is the Euler dilogarithm, Li2(l) = tt'^/6. 

In all the distributions, it is only the Standard Model part of the cross-section which is 
infrared singular for x ^ 1. The anomalous part of the electromagnetic vertex is proportional 
to the photon momentum and therefore infrared safe. 

To remove the cut, it is necessary to calculate the soft photon part of the bremsstrahlung 
contribution. Furthermore, the photonic vertex correction to r pair production has to be 
added in order to cancel the infrared singularity. The sum of soft and virtual corrections is 
known (see, e.g. []T3|) to read 



X 



C{v^ + a^)l2{ln^ 



ln(l - x) + 



4 



+ 2Li2(l)- 



(21) 



In sum, the 0{a^) part of the total cross-section for r^r (7) production is given by 



3 1 1 s \ 3 s 2 
3 — m — TT H -fl 

4 2 1 ml 64 ml ^ 



1 



--a-r 
4 



192 m 



2 T 



(22) 
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The terms proportional to may be found aheady in and the Standard Model contribution 
is also well-known. 

The existence of a non-vanishing electric dipole moment dr would lead to minor modifica- 
tions of the formulae. In the matrix element one has to replace 

dr = eF3(0). (24) 

There appear no interferences of d^ with the electromagnetic vertex or with a,- in the cross- 
sections. The contributions of the electric dipole moment may be obtained by the following 
replacement: 

4-4 + 4- (25) 

Since the interferences depend only on Or, a combined fit to both a-r and d-r becomes possible. 



3 Numerical Results 

The numerical calculations have been done with sin = 0.2320. 

Figure ^ compares the influence of the terms proportional to Qr and a^. An anomalous 
cross-section contribution of less than about 1.5 pbarn leads to a limit of about |aT-| < 5% from 
the a1 terms alone. The same measurement yields from the complete anomalous contribution 
a limit of about —5.6% < < 4.4%. Similarly, a limit of la,-! < 3% transforms into about 
—3.6% < a^- < 2.4%. Evidently, at \ar\ ~ 2% and below this value linear terms in a^- are 
indispensible for a correct interpretation of data. 

Figure |^ shows the photon energy distribution of the anomalous terms. The influence of 
the linear terms is more pronounced than in the integrated cross-section. Due to the (1 — x') 
dependence, the anomalous terms are proportional to but get damped at high photon 
energy due to the overall factor (3' . 

Here we have to discuss an important cut of the analysis. It is impossible to identify 
photons when emitted parallel to the radiating r. Thus, a cone has to be cut around the 
production angles 2 of the r leptons in the cms. These angles may be related by expressing 
the scalar products pi/c, p^k in two different ways 



Pi /2k = pi/2 



k^- 



1=1 



Pl/2t 



Pl/2 



1 -/3l/2 cos6'*/2 



(26) 



The momentum components are taken from ( [A.14|) . The velocities of r"*"/ in the cms are 



A/2 



4m, 



r W/2)^ 



V s I [{I +x')t {I ~x') 13' cos ef j' 
and the angles between r+/~ and the photon in the cms: 



(27) 



cos 6*^ 



/2 



It 13' cose 



P1/2I l[{l+x')Til-x')P'cos9]_ 



:28) 
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Figure 2: The contribution of the anomalous terms to (Thardi^)- 



In the appendix, (^) is inverted explicitly. 

In order to take advantage from this, one has to derive the double-differential distribution 
in s' and cos 6*: 



a 



dx'd cos di 



(3'C{ (1 - x') 
1 



2(1 - x' 



1 



+ 



X 



(1 



s 



6a' 



X r 



16 



/\2 



Ax'{v' + a") + (1 - cos" ^i) (1 - x) 



m 



{v' + a') 



(29) 



The next integration may be performed with account of the cuts on either numerically 
or, after the above preparations, also analytically. 

The result is shown in figure The terms are nearly not influenced, while linear terms 
are reduced drastically. This is due to their peaking behavior in forward direction; a behavior 
which is even more pronounced for the Standard Model QED contributions. Nevertheless, the 
cuts do not change the conclusions of our study. 

Finally, figure ^ shows the percentage of anomalous cross-section contributions compared 
to the Standard Model prediction for |ar| = 5%. If the minimal photon energy is larger, the 
weight of them rises (but the number of events goes down). For negative a,-, the rate is only 
60% of that for positive a^. 
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Figure 3: In the upper figure the two (solid and dashed) lines in the middle are contributions 
of the anomalous terms quadratic in for \ar\ = 5%. The two (solid and dotted) lines above 
are the net anomalous contribution for Gt = +5% and the two (solid and dotted) lines below 
for Gr — —5%. In the lower figure we plot the same curves for |a,-| = 1%. 
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E.,„ (GeV) 



Figure 4: The relative contribution of anomalous terms compared to the Standard Model pre- 
diction, crfi2rTi^)/^hard{^)> for fixcd , as a function of the cut on the minimal photon energy 
E ■ 

To summarize, based on an analytical approach, we have shown that linear contributions 
in ar are of numerical relevance if the limits on the anomalous magnetic moment of the r are 
better than a few per cent. Thus, they have to be taken into account with the present LEP 1 
data. A combined fit to both and becomes possible. 
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A Phase Space and Integrals 

Here we will present the parameterizations of phase space and four-momenta. We take the 
ultra-relativistic limit in the electron mass while m,- is taken into account exactly if not stated 
otherwise. 

The phase space depends on five independent variables besides s: The spherical angle of 
the photon in the cms, the spherical angle of the in the rest system of the r-pair, and 
x' = s'/s, s' being the invariant mass of the r-pair. 

With these variables, the phase space for the three particle final state may be expressed 
as a sequence of two two-particle phase spaces: 

^T- ^^Pl ^^^^2 d''^^- .4/. , . .X 
di = 77- -t; K- rTT-r^O (Ki + Ko — Pi — P2 — k) 

(27r)32p? (27r)'^2p^ (27r)32A;0 ^ 2 i^i 1^2 j 



dr^dr^dx', (A.l) 



with 



and 



J\{x',ml/s,ml/s) j3' 
dr^- = — dv5idcos6'i = -— d99idcos6'i, (A. 2) 

dFc = ^ -dip^AcosO^, (A.3) 



x' = ^^i±^, (A.4) 



- Jl-'^. (A.5) 
V x's 

The A function is defined as usually: 

A(a, b, c) ^ + + - 2ab - 2ac - 2bc. (A.6) 
The kinematical limits are: 

0< ^Pu^Pj <27r, (A.7) 

0< 9i,9^ <TT, (A.8) 

— ^ < x' < 1. (A.9) 
The limits on x' are related to those on the photon energy: 

0<E^< E^^ = ^1 - . (A.IO) 
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A dedicated construction of the four-momenta is crucial for the subsequent steps of the 
calculation. We determine all momenta in the cms. and begin with an ansatz for the photon 
momentum: 

g = ^(l-x',0,0,l-,T'). (A.ll) 

Then, we introduce the photon production angle 6^ as the angle between the three-momenta 
of photon and electron in the cms: 

ki = ^(l,sine^,0,-cose^) (A.12) 

k2 = ^ (1, - sin 0^,0, cos e^). (A. 13) 

All the four-momenta are independent of ip^. 

Finally, we have some freedom to choose the momenta of the r leptons in the rest system 
of the r-pair. Following the arguments given in Appendix B of |]16|, we get after a Lorentz 
transformation into the cms: 

Pi/2 = ^^=Fa(l - x') + 1 + x', ±bcosipi, ±bsm(pi, ±a(l + x') - 1 + x' j , (A.14) 

a = /5'cos^i, (A. 15) 

b = 2v^/3'sin0i. (A.16) 
The denominators of the r propagators are: 

U = 2p,q = |(l-x')(l-/3'cos^i), (A.17) 

t^ = 2p2q = |(l-x')(l + /?'cos^i). (A.18) 

The first integrations are over (pi and cos^^ and are simply polynomial. The only non- 
polynomial angular integrations are over 6i: 

1 1 2 / 

- / dcos^i^ = — ^-TT, (A.19) 



-1 

1 



2 7/^°^^^^ = /3'.(i-x') ^"r^"/5'.(i-x') ^""^^ ^^-^'^ 

The approximation in ( |A.2(j| ) is used only for the integration over x' where it does not lead 
to substantial inaccuracies. For the integrations with angular cut we have to invert (pSl) . 
Equation ( A.14| ) may be expressed explicitly in terms of cos6'^^2- 



1 

A/2 



Pi/2 = \Pi/2\ I ^5—, sin 91/2 cos ipl/2, sin 6^1/2 sin v^I/a, cos 6*1/2 > (A.21) 
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with 

Pl/2 



cos^r/2 = 7^- (A.22) 

' \Vl/2\ 



From pIi2 and IP1/2I = Jp^\/2 ~^ ~^ P^i/2 ^oi^ ( |A.14| ), we find immediately 



±Acos9i-B s 

cosg*/, = ^ , (A.23) 

' V^(l - cos2 9i) + A'^ cos2 9i T ^AB cos ^1 + 5^ 

with 

^ = /?'(! + x'), (A.24) 

5 = (A.25) 

D = Ax'P'\ (A.26) 

E = A^ - {A^ - D) cos^ 9*y2 (A.27) 
The inversion of (|A.23| ) leads to a quadratic equation for cos6'i/2 with two solutions: 



AB{1 - cos2 91) ± JA^B^{1 - cos2 91^ - E{B^ - {B^ + D) cos2 ^t) 
cos 6*1 = — (A. 28) 

In the limit x' — > 1 (no boost) the corresponding angles in both systems must be equal. 
Therefore, the sign in front of the square root has to agree with the sign of cos6'*y2. 

The cuts on cos 6** lead to a reduction of the range of integration over cos 6*1. With ( |A.28| ) 



the new integration limits are known as functions of cos^*. Since one has to apply cuts 
on both T leptons, the integration limits remain symmetric and the following modifications 
of (|A:T9|) -( [OOD arise: 



1 /dcos^i^ = 0(^)^0, (A.29) 
z J s 

— c 

1 r 1 1 1 + c 

-/dcos^i- = ^In-^. (A.30) 

2 7 t± (1 — x'js 1 — c 

— c 

The only potentially non-trivial integrals are those over x'. It is here where the limits 
/5' — 1 and m^/x' — lead to simplifications. Aside from integrals over (x')"", (x')" Inx', there 
are two integrals which require an infrared cutoff x: 

ln(x-3/m;) ^^, = ln^ln-^-Li2(l) + Li2(l-x), (A.31) 
1 — x' 1 — X 

^dx' = In^^. (A.32) 
1 - x' 1 - X ^ ' 
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